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Abstract 

The  structure  of  a  plane,  steady  hydromagnetlc  shock  Is 
studied  from  the  point  of  view  of  one -fluid  continuum  theory. 
Viscosity  and  heat  conductivity  are  neglected;  the  electrical 
conductivity  may  have  any  value.   The  purpose  is  to  derive 
expressions  for  the  shock  profile  and  thickness  in  terms  of 
the  parameters  that  enter  into  the  equations  of  motion  of  the 
fluid. 
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NYO-7693 


ON  THE  STRUCTURE  OF  A  STEADY  HYDROMAGNETIG  SHOCK 


One -Fluid  Theory 


Marian  H.  Rose 


I.   Introduction 

The  properties  of  hydromagnetlc  shocks,  considered  as  sharp 
discontinuities,  are  now  well  known  through  the  work  of  de  Hoffmann 
and  Teller  ,  Heifer  ,  Lust   and  Frledrlchs^.   However,  due  to  the 
relative  novelty  of  the  subject,  the  problem  of  the  Internal  struc- 
ture of  these  shocks  has  not  yet  received  much  attention.   Recently, 
Marshall  has  treated  the  case  of  a  plane,  steady  hydromagnetlc 
shock  with  constant  viscosity  and  heat  conductivity  and  has  obtained 
profiles  for  very  high  and  very  low  electrical  conductivities »   He 
shows  that  among  other  parameters  the  thickness  is  proportional  to 

1.  De  Hoffmann  and  Teller,  Magneto-Hydrodynaralc  Shocks,  Physical 
Review,  80,  692-703,  1950. 

2.  Heifer,  L. ,  Magneto -Hydro dynamic  Shock  Waves,  Astrophysical 
Journal,  llj,  177-199,  1953. 

3-   Lust,  R.,  Magneto-Hydrodynaraische  Stossweller  in  einem  Plasma 
unendllcher  Leltf ahigkeit » 

Zeltschrlft  fur  Naturforschiing  Band  8a,  277-281^.,  1953. 
Statlonare  Magneto-Hydrodynamlsche  Stosswellen  belleblger  Starke » 
Zeltschrlft  fur  Naturforschiing  Band  10a,  125,  1955. 

l\..      Friedrlchs,  K.O.,  Non-Linear  Wave  Motion  in  Magneto-Hydrodynamics, 
Mimeographed  Notes  -  Los  Alamos,  1955» 

5.   Marshall,  W. ,  The  Structure  of  Magneto-Hydrodynamlc  Shock  Waves  - 
Proceedings  of  the  Royal  Society,  A,  233,  367-376,  1956. 
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the  mean  free  path.   Essentially  the  same  results  have  been  obtained 
by  Sen  for  the  case  of  infinite  conductivity. 

In  this  paper,  we  shall  consider  a  steady  shock  in  a  medium 
where  the  state  variables  may  change  by  significant  amounts  over 
the  distance  of  a  mean  free  path,  i.e.,  the  mean  free  path  is 
assumed  to  be  very  large.   Consequently  heat  conductivity  and  vis- 
cosity lose  their  meaning.   It  is  also  assumed  that  the  fluid  is 
completely  ionized.   The  simplest  case,  the  one  that  will  be 
treated  here,  is  that  of  a  plane  shock  propagating  in  a  direction 
perpendicular  to  the  applied  magnetic  field.   Obviously,  a  more 
satisfactory  treatment  of  hydromagnetic  shocks  would  be  based  on 
a  two-fluid  theory  which  includes  dissipative  effects.   However, 
even  the  simple  picture  given  here  enables  one  to  gain  some  in- 
sight into  the  behavior  of  these  shocks,  and  interesting  proper- 
ties emerge. 

The  author  wishes  to  thank  Professors  Grad,  Kolodner  and 
Berkowitz  for  helpful  discussions  and  criticisms. 


6.  Sen,  Hari  K.,  Structure  of  a  Magneto -Hydro dynamic  Shock  Wave 
in  a  Plasma  of  Infinite  Conductivity.  Physical  Review,  102, 
5,  1956. 
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II.   Basic  Equations 

We  consider  the  case  of  a  completely  Ionized  fluid  traveling 
with  velocity  u  with  respect  to  the  shock  in  the  positive   x 
direction.   The  magnetic  field  is  in  the   z  direction.   All  the 
state  variables,  riamely  the  density   p^  the  pressure   p,  the 
fluid  velocity  u,  the  magnetic  field  B,  the  current   J  and  the 

electric  field  E  depend  on  x 


/h 


I 


°^^'  alone  (Figure  1).   Viscosity  and 

heat  conductivity  are  neglected 
(see  Section  I);  the  permeability 
[I     and  the  electrical  conductivity 
u(x)       cf    are  ass;imed  constant o   The 


->  X 


FIG.  1  motion  of  the  fluid  is  governed 

by  the  equations  of  conservation  of  mass,  momentiim  and  energy 
together  with  Maxwell's  equations  and  Ohm's  lawo   To  complete 
this  system  of  equations,  it  is  necessary  to  add  an  equation  of 
state  for  the  fluid  and  to  make  some  assumption  about  its  internal 
energy.   Here  we  shall  assuune  that  the  fluid  is  a  perfect  gas  of 
monatomic  particles.   Thus,  for  a  steady  state,  we  obtain  the 
following  set  of  equations 

(1)  Conservation  of  mass       pu  =  M 

(2)  Conservation  of  momentum  \7  ( P'ti  +  p)  =  JxB 

d       ou^ 

(3)  Conservation  of  energy  •jr[u(e    +     2~  "*"  P^    =  E»J 


-  6 


where   M  Is  a  constant  of  integration  and  e   is  the  Internal 
energy  per  iinlt  masso   The  symbol  \7   denotes  the   x   component 
of  the  gradient. 

Equation  of  state:   e  =  "IrT  where   T   is  the  absolute  temper- 
ature and  R   is  the  gas  constant o 
Maxwell's  equations: 
(l^.)         curl  E  =  0 

(5)  curl  ^  =  J 

(6)  div  B  =  0 

(7)  Ohm's  law:   J  =  o'(E  +  uxB) 

using  rationalized  M,KoSo  units.   We  note  that  the  Hall  term  in 
Ohm's  law  has  been  dropped  since,  for  this  geometry,  it  is  easily 
shown  to  vanish  identically  from  the  equations  into  which  J  enters,^ 
Combining  (i|),  (5)  and  (7)  and  then  integrating,  we  find 

where  /\  is  an  integration  constant.   Equations  (2)  and  (5)  yield, 
after  integration; 

2 
(2')  pu^+p+|j2  =  iT 

TT  being  an  integration  constant o   Lastly,  using  equation  (8)  in 
equation  (2)  for  the  conservation  of  energy  leads  to 

(3-)  ti(|  pu2  +1  p  +|-)  =1 

8.      Grad,   H.,    Notes   on  Magneto-Hydrodynamics   IV,    NY0-6i|86,    New  York 
University,    1956o 
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after  Integration,   b   being  another  integration  constant.   We 
note  that  in  this  last  equation  we  have  neglected  a  term  of  the 
form  ~l^'T~   '      This  approximation  is  certainly  justified  for 
weak  shocks  or  for  high  conductivities. 

Equations  (1),  (2'),  (3')  and  (b)  constitute  the  basic 
system  of  equations.   They  contain  four  Integration  constants 
which  are  the  parameters  of  the  fluid.   Fortunately,  this  niomber 
can  be  reduced  by  two  by  rewriting  the  equations  dlraensionlessly. 

For  this  purpose,  we  introduce  the  following  set  of  dimen- 
slonless  variables  and  parameters; 

,2 


p 

2^LTr' 

Mu 

TT 
P    =       p    P    , 

I 

_  ML 

/\  -       ^"^      A 

X    = 

M.cJ'TT             -   _  M^e 

L\        / — 3A 

M     ^'      °           ^2 

The  equations  become 

(9)  pu  =  1 

(10)  u  +  5  +  B^  =  1 

(11)  u(u  +  5p  +  i4i^)  =  £ 

(12)  i=Bu-A 

From  here  on,  remembering  that  all  quantities  are  dlmensionless, 
we  shall  omit  the  bar  above  the  various  quantities  since  there 
is  no  longer  any  possibility  of  confusion. 

There  now  remain  only  two  parameters,   £   and  [\    .      The 
parameter  t      occurs  even  when  no  magnetic  field  is  present  and 
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Is  closely  connected  to  the  shock  strength  parameter   X  of 
ordinary  fluid  dynamics.   For  instance,  one  possible  definition 
of  X  is^: 

X  =  s/2.S   -  16^ 

On  the  other  hand,  /\  is  an  entirely  new  parameter  which 
depends  on  the  amount  of  magnetic  field  present.   Since  there 
are  two  parameters,  two  auxiliary  conditions  should  be  given  in 
order  to  specify  the  problem  completely.   For  example,  the  values 
of  u  and  B  might  be  assigned  in  front  of  the  shock*   These 
two  conditions,  together  with  the  equations,  are  sufficient  to 
determine  the  values  of  the  remaining  variables  and  parameters, 
both  in  front  of  and  behind  the  shock.   It  is  now  a  simple  matter 
to  reduce  the  system  of  differential  equations  to  a  single  dif- 
ferential equation  for  u  which  may  be  integrated  numerically. 
However,  before  having  recourse  to  numerical  methods,  we  shall 
first  describe,  in  the  next  two  sections,  some  interesting 
properties  of  the  shock. 

In  the  last  sections,  a  method  for  obtaining  weak  shock  pro- 
files will  be  outlined.   One  of  the  most  important  aspects  of 
the  shock,  namely  its  thickness,  will  also  be  discussed  and  its 
dependence  on  the  various  parameters  will  be  shown  explicitly. 


9.   Grad,  H. ,  The  Profile  of  a  Steady  Shock  Wave.   Communications 
on  Pure  and  Applied  Mathematics,  Vol,  V,  257-300,  1952. 
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III.   Discussion  of  Equations 

Since   p  and   p  are  given  in  terms  of  u  and  B  by 
integrated  relations,  it  is  possible  to  express  the  equations 
in  terms  of  u  and  B  only.   The  differential  equation  (12) 
for  B   does  not  yield  any  information  about  the  domain  of 
existence  of  the  shock.   This  information,  however,  may  be 
obtained  from  the  differential  equation  for  u.   From  equations 
(9),  (10),  (11)  and  (12),  we  obtain 


(13) 


(11^.) 


^ 


r   dn  _ 

2Bu 

2(Bu   -  A) 

dx 

I     -  i,u2 

>       B^  = 

5  - 

ifu  -^ 

Of  all  the  possible  solutions  of  this  system,  the  one  which 
represents  a  shock  must  satisfy  the  following  conditions S 
u  must  be  a  single-valued  function  of  x  and  must  approach 
a  constant  value  as  x  approaches  infinity.   Letting  -r—  =   0 
entails  Bu  =  /\,      Substituting  for  B   in  terms  of  u  in  (13) 
leads  to  the  following  cubic  for  u; 


(15) 


"^-!"^^f-^  =  o 


Since  /\      >   0,   and  since  only  positive  roots  of  (1^))  are 
relevant,  the  third  extraneous  root  must  be  negative.   The 
two  roots  of  interest,   u   and  u  ,   will  then  represent  the 
values  of  u  at  infinity,  in  front  of  and  behind  the  shock. 
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These  statements  are  illustrated  in  figure  2.   The  roots,   u^ 
and   u  ,   lie  at  the  Intersection  of   Bu  =  /\     with  the  curve 


B'^  =  5  -  1+u 


u 


B. 


^PURE  HYDRODYNAMIC  ROOTS 


FIG     2 


-  -  corresponds  to  By  =A 
-  corresponds  to 


B=5-4u-—  or   Gcurve 


A  zero  strength  shock  occurs  when  the  two  curves  beccrtie 
tangent  to  each  other.   Then  u   =  u  =  \i    (£).   We  note 
that  the  points  of  intersection  of  the  "£  curve"  with  the 
u   axis  correspond  to  the  roots  for  the  pure  hydro dynamic 
case  and  that  the  roots  for  the  hydromagnetic  case  always 
lie  between  these  two. 
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It  is  immediately  apparent  that,  for  a  given   C,  Bu  -  /\     <   0 
where  /\   is  the  value  of  /\_  for  a  zero  strength  shock.   Hence, 
the  ni-unerator  of  -r-  is  always  negative  or  zero.   On  the  other 
hand,  the  denominator  changes  sign  at  u  =  /j-  =  u  (^)   where 
■^  becomes  infinite.   Consequently,  for  u  </t-  ,   a  shock  is  no 
longer  possible. 


10, 


The  value   ^  ~  /t7  also  corresponds  to  the  maximum  of  B   = 

5  -  I4-U  -  —  .   Hence,  no  amplification  of  the  magnetic  field 

within  the  shock  is  possible.  The  opposite  statement  had  pre- 
viously been  made  (cf.  Rose  and  Grad,  Bulletin  of  the  American 
Physical  Society,  Series  II,  Vol.  I,  No.  1,  1956). 
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IV.   Ranp;e  of  Parameters 

It  is  apparent  from  the  preceding  discussion  that  shocks 
are  obtained  only  within  certain  ranges  of  values  of  the  parameters, 
Before  determining  these  ranges,  we  shall  introduce  a  new 
parameter  -y\.    ,      related  to  the  shock  strength.   For  reasons  of 
convenience,  we  shall  consider  ^  and   _A_  to  be  the  fundamental 
parameters  of  which  /\^     is  a  known  functiono 

In  the  pure  hydrodynamlc  case,    the  values  of  u,   far 
from  the  shock,  are  the  roots  of  a  second  degree  equation,  namely, 

u+  -^^   (5  +/25  -  16^)   - 

An  obvious  choice  for  the  shock  strength  parameter  is  then 

X  =  ^25  -  16  E)   or  some  convenient  fxinction  thereof,  since  for 

zero  shock  strength,   u  =  u   (u  =  •§  in  this  case)., 

In  an  analogous  way,  the  shock  strength  parameter  for  the 
hydromagnetic  case  is  defined  from 


u  =  u 
o 


+  f(A,  £  ) 

where   f (vL,  6  )   may  be  determined  when  -/V  and   £-  are  given. 
Some  information  about  the  range  of  the  parameters  may  be 
obtained  from  the  conditions  for  a  zero  strength  shocko   A 
necessary  condition  for  such  a  shock  is  that  'T"  =  0   in  equation 

11.   Cf„  reference  9. 
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(13),  leading  to  equation  (l5)  for  u.   Since,  furthermore,  we 
seek  two  roots  of  this  equation  which  are  equal,  positive  and 
real,  the  third  root  must  also  be  real.   Hence  (15)  may  be  ex- 
pressed in  the  form 

(u  -  u^)^  (u  +  u^)  =  0 

leading  to  the  following  relations  between  the  coefficients: 


(16)  -  t  =  ^o  -  2% 


(18)  4^  =  u'u^ 


I|.      o  o 


Prom  (16)  and  (17) 
(19) 


2 


Prom  (18),  we  know  that  u'  >  0.   Substituting  for  u   in  (18) 


shows  that  the   +   sign  in  the  expression  for  u   must  be  taken 

to  obtain  a  positive  value  of  u'.   Hence 

o 


12 
From  (16)  we  obtain 


u:  =  -  |p  +  y25/36  -  a/3 


(20)     A^  =  C-  I2  ■"  >/2^/36   -    £  /3  )    (5/6   +725/36   -    b/3)-  =  Ao 

Both     u        and     /\        depend  on      £     alone.      Since      £     may  vary 
between  limits,    which  will   presently  be    determined,    we   see    that 
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there,  exists  a  one-parameter  family  of  zero  strength  shocks.   As 
we  should  expect,  the  zero  strength  values  for  pure  hydro dynamical 
shocks  are  contained  within  this  family,  namely,  for   ^  ~  16 
we  obtain  u  =  ^  and  A  =0.   The  "£  curve"  then  loses  its 

o    O         *— 'O 

structure,  the  points  of  intersection  u.^,  and  u_  with  the   u 
axis  coalescing  at  the  point  u  =  ■&.   In  this  case,  therefore, 
only  a  zero  strength  hydromagnetic  shock  is  possible.   From  (19) 
it  is  obvious  that  the  upper  limit  of   C   is  j^.   The  lower 
limit  is  obtained  by  setting  p  =  0  in  the  equations  of  motion 

far  from  the  shock  (p  may  not  be  negative).   Thus  equation  (10) 

2  2 

becomes   B  =  1  -  u  and  substituting  for  B   in  (II4.)  yields 

c^  =  ij-u  -  3u  .   The  minimum  value  of  c   may  be  obtained  by 

differentiating  this  expression  with  respect  to  u.   The  result 

is   £  =  ■^.   The  parameter  C  must,  therefore,  lie  within  the 

range 

On  the  other  hand,  y\_  may  vary  from  0  to  -A-  where 

<'  </  max 

y\^  is  the  maximum  possible  shock  strength  for  a  given   c  . 

This  occurs  when  u   is  chosen  to  coincide  with  u  ,  i.e., 

A  =  u   -  u  .   Hence  -^   lies  within  the  range 

V  \_  max    o    c  ^ 


0  <  A  <  V^  -^/2^/36  -  t/3     . /E     . 


In  figure  3i  the  permissible  values  of  j\-    and   B  lie  within 

the  shaded  area.   We  note  that  -/>-     is  a  monotonlcally 

max  '' 
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decreasing  function  of   c  , 
the  maximum  value  oc curing 
at  £  =  ^  for  which  -A 
=  O.Qb9.      This  value  of  -A- 
corresponds,  therefore,  to 
the  strongest  possible  shock, 
The  range  of  /\     may  be 
determined  in  an  entirely 
similar  way.   It  is 


Ao  <  A  <  Zio 

where  A    ^^  defined  by  equation  (20)  and  Ac  ~  ^c^c 

For  purposes  of  future  reference,  the  expressions  for 
^(X,£(u  ))  and  f(A,£(u^))  will  be  given  here.  Prom 
/\   =  u  B  ,   and,  using  equations  (17)  and  (It-)  to  express 


£.   in  terms  of   u  ,   we  obtain 


(21)  A(A,£(u^))  =  (u^-yV)[5-^(%-A)   - 


^^.^1  -  3u^) 

-o  -A 


o  2 


1/2 


] 


In  order  to  determine   f(A,  £(u^)),   we  allow  A   ^o  vary 
about  its  zero  strength  value  in  equation  US)    while  keeping 
(S   fixed  {  £    =    £.     ).   The  solutions  will  also  vary  about 
their  zero  strength  value  and,  therefore,  must  satisfy 


[u-  (u^-A)Hu-  (u^+  f(A,  e  ))][u+  u^  +  TVI 


0 
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Co3,-aarison  of   coefficients   leais   to 


..,  -  /  -  i' 


(22)      fCA,  f  Cu„))    =   -  ^^^ ^U   -   /l   -^       ^-^~ 


3-a-^-  t 


J 
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V.   Shock  Profile 


We  shall  here  describe  a  method  for  obtaining  weak  shock 
profiles  from  which  values  of  the  thickness  may  be  deduced. 
The  method  consists  in  assuitilng  the  following  form  for  -r—  ;  namely. 


(23) 


1^  -  I  (u,A,  £)(u  -  u^)(u  -  u_) 


and  then  expanding  K  (u,A  ,  6)   in  a  power  series  about  its 
value  at  yV  =0,   Thus,   K(u,/\,,  £  )   may  be  expressed  as 

(21^.)   K(u,A,  £)-K    (u,6)+AK    (u,  £)+Ak    (u,e)  +  ... 

The  coefficients  of  j\.    are  determined  by  comparison  with  equation 
(13).   The  function  K(u,A,£)   has  some  physical  significance 
since  it  can  be  shown  to  be  related  to  the  shock  thickness  f  . 
Although  TT  may  be  defined  in  a  variety  of  ways,  for  our  pur- 
poses, the  following  definition  is  meaningful} 


US) 


V 


u  -  u , 


1—1 
'dx' 


max 


where 


,du 

dx 


is  the  maximum  slope  in  the  profile  (cf.  figure  l\.)  „ 


max 


Assuming  ^  constant,  X'    would  have  the  value  TT  =  ~ — ' — 


u  -  u 


FIG.  if 


.  _1L 
->  X 


K(u  -  u^) 


Hence,   K  would  be 
inversely  proportional 
to  the  thickness.   In  this 
case,  equation  (23)  inte- 
grates to 
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2x       -1  2^-'^-  "^  ^+^ 

(26)  -  ^  =  tanh  "^  — 

'^  u   -  u 

which  satisfies  the  appropriate  conditions  at  x  =  +  gd  . 

Before  proceeding  with  the  calculations,  it  is  necessary 
to  make  a  change  of  variable,  the  reason  being  that  the  pro- 
file degenerates  to  a  straight  line  for  j\-   =  0  ;   u  therefore 
cannot  be  analytic  in  JV     in  this  neighborhoodo   By  intro- 
ducing the  proper  variable,   w,   the  profile  retains  its 
structure  even  for  zero  strength  shocks.   We  define   w  as 
follows; 

(27)  u  =  u^  +  TLw 

hence,   -1  <  w  <  ^      \'  where  ^''}y       ^   — >  1   when  A.—>   0. 

We  also  introduce  a  new  coordinate  y  which  may  be  obtained 
from  X  by  the  transformation  y  =  ^Xo  Equation  (23)  then 
becomes 


(28) 


dW  _  T7 


=  K(w,A,  <i  )(w  +  l)(w  =.  ^SA^U.)    , 


Similarly,  equation  (13)  transforms  into 


dw 


2(u^^Aw)2  /5  -  -^-^   -  i,(u^-H  Aw) 

_  V        O 


A  [£  -4(u^+ Aw)''] 


(29) 


(u  +  Aw) 


7^- vfe"''^"°"~^"^ 

_(o)  _(1)       _(n) 

K    ,K    ,«„,,K    ,000   are  then  determined  by  comparison 
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of  equations  ( 2ti )  and  (29),  the  latter  first  having  been 
expanded  in  a  Taylor  series  about  u  .   Since  the  main 
purpose  of  this  paper  is  to  determine  the  general  behavior 
of  the  shock  thickness  rather  than  to  compute  the  profile  to 
higher  order  terms,  we  shall  limit  ourselves  to  the  first 
term  in  the  expansion.   The  result  is 

(30)  K    =  2_^  U(|  -  u  )  +r^(7  -  pi 

2(5-8u^)[  2    o    i;      2u 

a  quantity  independent  of  w.   In  this  approximation,   w  as 
a  function  of  x  is  simply 

(31)  w  = ^ tanh  j^^ K   y  t    VL  '  '  ^ 

hence,         / 

(32)  u  =  u^+A ^^^k_tanh  ^-j^     ^^  ^^J^  '  ^.  Mx  +    VL 

The  profile  of   B   is  obtained  by  substituting  for   u   its 
above  value  in  equation  (li|).   In  figures  5,  6,  7  and  8,  pro- 
files of  u   and   B   are  drawn  for  various  values  of  VL  and   c« 
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VI.   Shock  Thickness 

Using  equation  (25),  we  find  that  the  value  of  the  shock 
thickness,  to  a  first  approximation,  is 

(33)  r  =       ^ ,  . 

(The  above  value  for  ^     has  been  given  the  dimensions  of  a 

M 
length  through  multiplication  by  ~"^)«   Replacing  M  and  rr 

by  their  values  from  equations  (1)  and  (2')  and  assuming,  as 

3/2 
a  rough  approximation,  that  C' oaT-^'     ,   we  conclude  that 

1)  7^  oc  — n   keeping  p,   T  and  u  fixed 

2)  r  00  T-^/2  ^^^  P  >  |i 

p 
t>  oc  T"-''    for  P  ^'^  "577  keeping  p,  B  and  u  fixed  , 

It  is  worthwhile  comparing  the  shock  thickness  with  the 
Larmor  radius  L  which  might  plausibly  replace  the  mean  free 
path  as  a  limiting  thickness  for  the  shock  when  the  mean  free 
path  becomes  very  large  compared  to  the  dimensions  of  the 
apparatus.   Now,   ^  =  Tg  where   m  Is  the  mass  of  a 
particle,  and  v,   is  its  average  velocity  (v  oc\/T  where   T 
is  the  temperature  of  the  gas).   Hence, 

^  k  pueB 

-   -/^    A         tiCi'mv  (Mu  +  p  +  1^  ) 

Dividing  niimerator  and  denominator  by      p     leads   to   the   expression 
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£  +1 


2 


■  o    in  the  denominator  which  correspoads  to  the  square 
of  one  of  the  three  possible  speeds  at  which  a  hydromagnetic 

shock  may  propagate.   These  have  been  termed  "slow",  "inter- 

12 
mediate"  and  "fast"  by  Priedrichs   .   For  o\xr   geometry,  the 

shock  can  be  shown  to  propagate  only  at  the  fast  speed, 

namely,   u„   ^  =  /  ■"  +  p-  •  «   Assiiming  that   u  is  approximately 

equal  to  u„   ,   we  obtain 


b2 
For  small  values  of  -r^ —  ,  ^  is  roughly  proportional  to  B 

for  a  given  T   (the  conductivity   0*  depends  on  T  but  is 

independent  of  B),   For  -x —  »  1,   t-  becomes  independent 

of  B.   In  this  case,  therefore,  V      and  L  are  proportional 

to  each  other  for  a  given  T.   In  other  worda  the  Larmor 

radius  apparently  represents  a  limiting  thickness  of  the 

shock  when  the  mean  free  path  is  infinite.   These  results  are 

siimmarized  in  figure  9  which  shows  how  t-  depends  on  B. 


12o   Cf.  reference  I4.. 
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VII »   Comparison  with  Marshall's  Results 

In  conclusion.  It  Is  Interesting  to  compare  our  results 
with  those  of  Marshall.   His  equations  Include  the  effects 
of  viscosity  and  heat  conductivity  and  are  therefore  appli- 
cable to  a  medium  where  the  gradients  are  small  compared  to 
the  mean  free  path.   On  the  other  hand,  the  case  considered 
here  Is  one  where  the  mean  free  path  Is  large  compared  to 
the  dimensions  of  the  region  under  observation.   Hence  the 
usual  concepts  of  viscosity  and  heat  conductivity  lose  their 
meaning.   In  spite  of  these  differences,  we  shall  show  that 
the  two  theories  do  Indeed  meet  at  a  certain  point.   For  this 
purpose,  it  will  be  helpful  to  recall  briefly  some  of  Marshall's 
results.   He  shows  that  the  nature  of  his  solution  depends 

on  two  dlmenslonless  parameters,  namely,   a  =  t^^     and 

3dc^ 
(3  =  £^^g.   where  \     and  -r^  are  the  coefficients  of  heat 

conductivity  and  viscosity  respectively,  and  c,   the  sound 

speed  of  the  undisturbed  gas.   Marshall  obtains  profiles  for 

(3  «  1  (cr  large)   and   p  »  a  (o"  small).   The  first  case 

cannot  apply  to  ours  since   p  »  1   when  t^  ' — >  0   for  any 

value  of  C   On  the  other  hand,  the  condition   (3  »  a  is 

always  met  for  any  6~    when  both  A  and  ^J,  approach  zero. 

In  this  case,  Marshall  shows  that  two  solutions  are  possible; 

for  one,  the  fluid  speed  behind  the  shock   is  greater  than 

the  local  sound  speed,  for  the  other,  it  is  less.   However, 

13.   Marshall  uses  a  coordinate  system  in  which  the  shock  is 
at  rest, 

-  2h  - 


it  is  easily  seen  that  only  the  direction  field  of  the  first 
solution  may  coincide  with  oxors.   The  shock  thickness,  in 
this  case,  is  proportional  to   p/ =  rr^  where  J^  -  \~     is 
of  the  order  of  magnitude  of  the  mean  free  path  in  the 
undisturbed  gas.   This  is  the  same  result  (except  for  pro- 
portionality factors)  that  was  derived  in  equation  (32). 

The  important  difference  between  the  two  cases  is  that 
Marshall  is  limited  to  wide  shocks  since   cr  is  small,  whereas, 
in  our  case  the  shock  may  become  narrow  since   0^  may  have 
any  value . 
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